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For thermoelectric power generation in a multi-terminal geometry, strong numerical evidence for a 
universal bound as a function of the magnetic-field induced asymmetry of the non-diagonal Onsager 
coefficients is presented. This bound implies, inter alia, that the power vanishes at least linearly 
when the maximal efficiency is approached. In particular, this result rules out that Carnot efficiency 
can be reached at finite power, which an analysis based on the second law only, would, in principle, 
allow. 
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I. INTRODUCTION 

Thermoelectric generators use a heat current between 
two reservoirs of temperatures Th and Tc < Th connected 
via a device to drive a particle current against a gradient 
in chemical potential, thus delivering useful power. As a 
fundamental consequence of the second law of thermody¬ 
namics, the efficiency 77 of such a heat engine is bounded 
by the Carnot value r]c = 1 — ^c/^h- This bound can 
indeed be reached if the device acts as a perfect energy 
filter ( 3 , however, only at the price of vanishing power 

ii- 

The relationship between thermoelectric power and ef¬ 
ficiency beyond this limiting case has been widely studied 
during the last years , see also 0 and [? ] for re¬ 
cent reviews. A particular landmark was achieved by 
Whitney [Hilll , who calculated the transmission func¬ 
tion that leads to optimal efficiency at given power and 
thereby identified upper bounds on the output power of 
quantum coherent devices. In all cases, however, finite 
power is only possible at an efficiency strictly smaller 
than ? 7 c- 

This antivalence resembles the power-efficiency- 
dilemma, which has been studied by Allahverdyan et al. 
on the basis of a quite general model for a reciprocating 
quantum heat engine [1^ . Their results confirm the gen¬ 
eral expectation that, under realistic conditions, a con¬ 
ventional cyclic heat engine can approach Carnot effi¬ 
ciency only in the quasi-static limit, within which its 
power output necessarily goes to zero. 

In strong contrast, Benenti et al. have recently pointed 
out that a magnetic field, which breaks the symmetry 
between the off-diagonal Onsager coefficients, might en¬ 
hance the efficiency of a thermoelectric heat engine sub¬ 
stantially such that rjc seems to be achievable even at 
finite power [T5|. This exciting suggestion follows from 
a straightforward analysis of the linear response regime 
using only the second law. However, so far, neither a spe¬ 
cific model for a thermoelectric generator delivering finite 
power at Carnot efficiency, nor a fundamental principle 
forbidding the existence of such a super device has been 
discovered. 

A promising platform for the investigation of both of 
these aspects is provided by the multi-terminal setup 


sketched in Fig. [H By operating this model as a thermo¬ 
electric heat engine, we previously derived bounds on its 
efficiency 77 that are stronger than the ones obtained by 
Benenti et al. 0 II 3 . These bounds, which follow from 
current conservation, depend only on the number n of 
terminals and, for any finite n, constrain 77 to be strictly 
smaller than 77 c whenever the off-diagonal Onsager coef¬ 
ficients are not identical thus preventing finite power at 
Carnot efficiency. 

In this paper, we move an essential step forward by 
investigating whether not only efficiency but also power 
can be bounded. Since power P carries a physical di¬ 
mension, such an endeavor can be expected to be harder, 
and, a priori, less universal, than searching for a bound on 
the dimensionless efficiency. Bounding power is arguably 
more relevant since a high efficiency is useless from a 
practical point of view if it comes with minuscule power. 
Here, we provide strong numerical evidence that P is in¬ 
deed subject to a new bound. This bound implies that 
the power output of the device vanishes whenever its ef¬ 
ficiency approaches the upper limit corresponding to the 
respective number of terminals. Finally, by extrapolating 
our findings to arbitrary n, we derive a conjecture for a 
universal bound on P, which rules out the option of finite 
power at Carnot efficiency even in the limit n ^ 00 . 

The paper is organized as follows. In section H, we give 
give a brief review of the essentials of the multi-terminal 
model and the bounds on its efficiency as a thermoelectric 
heat engine, which were derived in our previous work 
00 - Section HI contains our main results. We outline 
the numerical procedure used to substantiate our new 
constraint on the Onsager coefficients and show how this 
constraint can be used to bound the power of the multi¬ 
terminal thermoelectric generator. Finally, we conclude 
in section IV. 


II. MULTI-TERMINAL MODEL 
A. Setup 

We consider the model sketched in Fig. [H On the pri¬ 
mary level, n electronic reservoirs of respective tempera¬ 
ture Ti and chemical potential i^i (i = h, I,..., 77 , — 2, c). 
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FIG. 1. Sketch of a thermoelectric generator. Both, heat and 
particle current run from the hot reservoir on the left to the 
cold reservoir on the right thereby traversing the device D, 
which is represented by the dashed box. Within the multi¬ 
terminal setup, the device consists of a central scattering re¬ 
gion subject to a magnetic field B and n — 2 probe terminals. 
The figure corresponds to n = 4. 

are connected via perfect, one-dimensional leads to a cen¬ 
tral scattering region penetrated by the magnetic field B. 
The particle and heat currents flowing from the reservoir 
i to the scattering region are denoted by and J*, the 
corresponding affinities by 

= and (1) 

respectively, where the reference values /i = /ic and 
T = Tc have been chosen. In the linear response regime, 
to which we will stick throughout the paper, the phe¬ 
nomenological equations describing the coupled heat and 
particle transport in this system read 

j* = ^ Uj:f^ (2) 

jf=h,...,n —2 

with = (J;, J*)^ and hij e for 

i,jf = h, l,...,n — 2. Since, at this stage, we assume 
uncorrelated particles, which move coherently through 
the leads and the scattering region, the matrix blocks 


hij can be calculated explicitly using the multi-terminal 
Landauer formula [l8|, [l 9 | 

Uj=ldE f{E) {E-llf) 

with f{E) = cosh“^ {{E — ji)/{ 2 kBT)) /{AkBh)^ where h 
denotes Planck’s constant and kB Boltzmann’s constant. 
The energy and magnetic field dependent transmission 
coefficients Tij{E) encode the properties of the scatter¬ 
ing region and cover the effects of particle-particle in¬ 
teractions on the mean field level [^. Due to current 
conservation and time reversal symmetry, these dimen¬ 
sionless quantities must fulfill the sum rules [3 El 

TijiE)= = i ( 4 ) 

i=h,...,c jf=h,...,c 

and the symmetry relation 

Tij{E,B) = Tji{E,-B), ( 5 ) 

respectively. Throughout the paper, we notationally sup¬ 
press the dependence of any quantity on B, if there is no 
need to indicate it explicitly. For vanishing magnetic 
field, relation ([ 5 j) implies that the primary Onsager ma¬ 
trix showing up in (| 2 ]) must be symmetric. For a fixed 
B ^ 0 , however, (| 5 j) does not lead to any further con¬ 
straint. 

The full multi-terminal setup described above becomes 
a model for a thermoelectric heat engine by imposing the 
boundary conditions = 0 for i = 1 ,..., n — 2. 

Using these constraints to eliminate the corresponding 
affinities we obtain the reduced set of phenomeno¬ 

logical equations 



where the with (3 = p^q denote effective Onsager 
coefficients, which are given by 


f Lpp 

\Lqp 


Lpq \ 

Lqq J 


^Lhi,.. 

• : IL'hn —2^ 

/ Ln • 

Lln-2 \ 

-i 

/ Lih \ 



S 

1 

to 

■ • Ln-2 n-2/ 


\f^n—2h / 


( 7 ) 


with hin the blocks of primary Onsager coefficients de¬ 
fined in daj). 

The physical intuition behind this procedure is that 
only the hot and the cold terminal are considered as 
real, while the remaining ones act as probes, which do 
not exchange any net quantities with the physical reser¬ 
voirs. Originally proposed by Biittiker [^, such probe 
terminals have meanwhile become a well established 
method to phenomenologically model inelastic scattering 


events like electron-phonon interactions, see, for exam- 
ple, [3, El MM- Such processes would be hard to 
take into account on a microscopic level, but are, inter 
alia, crucial to obtain a non-symmetric Onsager matrix 

1313. 
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B. Bounds on efficiency 


on the Onsager coefficient Lqq, where 


The thermodynamic efficiency of a thermoelectric gen¬ 
erator in the linear response regime generally reads 


_ P _ _ ffi ^pql 

Jh- + ’ 


(8) 


where 


P = (Ai - = -T(ph)2 ^ ^ (g) 

is the supplied power, 

7 = ^pV-7 (10) 


denotes the ratio of the affinities and rjc = (Th — T)/T = 
TFq the linear response expression of the Carnot effi¬ 
ciency. Maximizing rj with respect to 7 yields the maxi¬ 
mum efficiency 


7max((^5 y) — 




^VyTT+1 

Here, we introduced the dimensionless parameters 

and , 

-^qp -^pp-^qq -^pq-^qp 


(11) 


( 12 ) 


which, due to the second law, have to obey the inequali¬ 
ties [iBl 

h{x) <^<0forx<0, 0<^< h{x) for x > 0 (13) 

with h{x) = 4x/(x — 1)^. For B = 0, x assumes the sym¬ 
metric value 1 as a consequence of time reversal symme¬ 
try and y reduces to the conventional figure of merit ZT 

In our previous work on the n-terminal setup [l7|, we 
showed that current conservation implies that m can 
be strengthened by replacing h{x) with 


hn{x) = 


Ax cos^(7r/n) 


{X - 1)2 • 

The resulting constraint leads to the bound 


7max((^) — ycX 


y/hn{x) + 1 — 1 
yjhn{x) + 1 + 1 


(14) 


(15) 


on the efficiency as a function of the asymmetry pa¬ 
rameter X, which, for any finite number of terminals, is 
stronger than the one required by the bare second law. 
However, no predictions on a putative maximal power 
can be drawn from the analysis presented so far. 


III. BOUNDS ON POWER 

The fundamental bound on power discussed in this pa¬ 
per follows from the constraint 


Lqq — Lqq/Nqq ^ 1 yjh^i^X^ 


(16) 




klT^q 


4h 


(17) 


is a dimensional factor with q defined in (|20D. We find 
this constraint, which constitutes our first main result, by 
randomly generating transmission matrices obeying the 
sum rules (|4|) and exploiting the matrix structure ([3j). 


A. Numerical procedure and results 

Our aim for this subsection is to numerically confirm 
the constraint m- Since it is convenient to work with 
dimensionless quantities, we first change the integration 
variable in © diS e = {E — y)/{kBT) and, second, define 
the rescaled matrix blocks 


u, ^ ly m yy {s„ - T.Ae)) 


(18) 


= eLijD, 


where 


T V 0 


(19) 


f{s) = cosh ^(^/ 2 ) and Tij{e) = TijisksT + y). The 
scaling parameters p and q are defined by 


P \ _ 


f 


de /(e) I ^2 


4 

47r^/3 


( 20 ) 


where u = y/ksT [^ . 

Next, to obtain a tractable parameterization of the ma¬ 
trix blocks m, we make use of the fact that the matrix 
T{e) with elements (¥( 5 ))^^ = Tij{s) is bistochastic, i.e., 
for any fixed 5 , its elements fulfill the same sum rules (|1]) 
as the %j{E). Consequently, by virtue of the Birkhoff- 
von Neumann theorem, there is a set of N permutation 
matrices and positive numbers \k (s) such that 


N 


N 


T(e) = E A/c( 5)P^ and ^A/c(5) = l. (21) 


k=l 


k=l 


Inserting this decomposition into m and formally car¬ 
rying out the integral yields 


N 


L.i=E 7 i- 7 ) 


( 22 ) 


k=l 


where the = (P^)ij are the matrix elements of the 
respective permutation matrix and 


( 23 ) 
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FIG. 2. Scatter plots of the rescaled Onsager coefficient Lqq 
as a function of the ratio yjhnix) for respectively 50000 ran¬ 
domly chosen n-terminal models. The dashed line is given by 
Lqq = 1 — y/hn{x). To generate these data, we used the set 
of permutations {(1)(2, 3), (1, 3)(2), (1, 3, 2)} for n = 3 (left 
column) and the set {(1)(2)(3, 4), (1, 4)(2, 3), (1,4, 2, 3)} for 
n = 4 (right column). While for all plots the sign of the Wk is 
chosen randomly, \wk\ is sampled from an increasingly large 
interval A C [0,1]. Specffically, we have from the hrst to the 
fourth line A = [1,1], A = [0.99,1], A = [0.5,1], A = [0,1]. 
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FIG. 3. Scatter plots of the rescaled Onsager coeffi¬ 
cient Lqq as a function of the ratio y/hn{x) for respec¬ 
tively 50000 randomly chosen n-terminal models. For 
all plots, the sign of the Wk is chosen by random and 

their absolute value has been set to = 1. The 

sets of permutations used to generate the shown data 

are {(1, 3)(2)(4, 5), (1)(2,4)(3, 5), (1, 4,3, 2, 5)} for n = 5, 

{(1)(2)(3,6)(4,5),(1,6,4,2,5,3),(1,6,4)(2,3,5)} for n = 6, 
{(1, 3)(2,4)(5)(6, 7), (1, 5,7)(2,6)(3,4), (1,7,4, 5, 2,6,3)} for 
n = 7 and {(1, 6)(2,8)(3,4, 7, 5), (1, 7,6, 3, 2, 8, 5)(4), (1,5,4) 
(2,3)(6, 8)(7)} for n = 8. 


is a positive semidefinite, symmetric matrix of dimension 
2. Due to these properties of the Mj,, there exist some 
numbers (Jk,ak G (0,oo), Wk G [—1,1], such that 


Mfe = Gk 


( ^ 

\Wkak 


Wkak\ 

4 )■ 


(24) 


Combining this expression for M/^ with the decomposi¬ 
tion ([22]) gives a parameterization of the dimensionless 
matrix blocks hij in terms of N permutation matrices 
of dimension n and 3N real parameters ak^ak^Wk- The 
latter are constrained by the two important sum rules 


N N 

k=l k=l 


(25) 


which follow directly by comparing (|M|) with (|23|) and 
using the sum rule m for the Xki^)- We note that, in 
principle, there is a third sum rule 


N 

(JkWkak = 

k=l 



= ln(16) + 4 In (cosh(z^/2)) — 2z^ tanh (z^/2). 


(26) 
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However, for simplicity, this constraint will not be ex¬ 
ploited within our analysis. 

We now proceed as follows. First, we chose a fixed 
set of N distinct, n-dimensional permutation matrices 
Ap = randomly pick a large number of pa¬ 
rameter sets Am = {cr/c, a/^, such that for any of 

these sets the sum rules (|25]) are fulfilled. Second, for 
any of the sets Am, we evaluate the matrix blocks hij 
according to and (I 22 | and subsequently calculate 
the rescaled effective Onsager matrix L using © with 
the hij replaced by hij. Third, for any individual of the 
thus obtained matrices L, we determine the values of the 
parameters x and y by inserting the elements of L into 
the definitions given in ([T 2 ]). This step is justified, since 
it is readily seen that the effective Onsager matrix L is 
connected to its rescaled counterpart via the transforma¬ 
tion L = and therefore the scaling factors con¬ 

tained in D cancel if the elements of L are plugged into 
m- Finally, we plot the rescaled Onsager coefficient 
Lqq, i.e., the lower right entry of the matrix L, against 
the ratio yjhnix). This choice of variables is quite rea¬ 
sonable, since, as we will show in the next subsection, 
the maximum power of the multi-terminal thermoelec¬ 
tric generator is proportional to Lqq = NqqLqq and we 
are especially interested in the behavior of this quantity 
as y approaches its bound hn{x). 

We begin with the minimal cases n = 3 and n = 4. 
Fig. [2] shows the results of the procedure outlined above 
for two representative sets Ap of respectively three per¬ 
mutation matrices and all Wk randomly chosen from an 
increasingly large interval. We observe that for any of 
these models the inequality m is respected and, for 
some of them, even saturated. Our numerical data fur¬ 
ther suggests that this bound is independent of the choice 
and number N of distinct permutation matrices in the 
set Hp, where N must be larger than two and at least 
one of the has to be non-symmetric to obtain a non- 
symmetric effective Onsager matrix. However, if N is 
increased, it becomes more and more improbable to find 
models, for which Lqq attains the bound (pT|) or comes 
even close to it. Likewise, we find that as the interval 
A, from which the Wk are drawn, is enlarged, the sharp 
boundary visible in the plots shown in Fig. [2] deterio¬ 
rates rapidly. Therefore, we are convinced that the data 
presented in the first row of Fig. [2] cover a representa¬ 
tive subset of the extreme points in the given parameter 
space with respect to the inequality m- 

For n > 4, in Fig. [3l we show representative data, 
which have been obtained for A = 3 permutation ma¬ 
trices P/c and \wk\ = 1. The bound ([T 6 |) holds for any 
of the randomly chosen models. However, in contrast to 
the cases n = 3 and n = 4, we were not able to achieve 
saturation. Even by increasing the number N beyond 3, 
this finding persists. Sampling the absolute values of the 
Wk from a finite interval A, like for 3 and 4 terminals, 
leads to a decay of the sharp boundary lines appearing 
in Fig. [3l It therefore remains an open question at this 
stage, whether or not, for n > 4, there are models that 



X 

FIG. 4. Bound (1281) on the rescaled maximum power output 
Pma^{x)/Po of the multi-terminal thermoelectric generator as 
a function of the asymmetry parameter x for n = 3, n = 4 
and n oc. 


saturate the bound m- 


B. Maximum power 

We will now demonstrate in two steps that the new 
bound (HH) significantly restricts the power supplied by a 
multi-terminal thermoelectric generator. First, we max¬ 
imize ([9]) with respect to 7 thus obtaining the maximum 
power 


Pmax(a:,y) = (27) 

Second, replacing Lqq by its upper bound ([T6|) . as indi¬ 
cated from now on by a hat, and optimizing the resulting 
expression with respect to y gives the upper bound [29| 


^(x) — Pinax{^,y (^)) — 


yjhn{x) + 1 — 1 
\/hn{x) + 1 + 1 


(28) 


where 


Po = T{T'^fN,jA (29) 

sets the scale and y*{x) = ^JhJ^x)~L^ — 1. The bound 
(|28)) . which is our second main result, is plotted in Fig. 01 
In the limit x Too, Pmax asymptotically reaches the 
value Pocos^( 7 r/n) and, irrespectively of n, Pmax(^) at¬ 
tains the global maximum Pq for x = 1. Moreover, Fig. |4] 
reveals that the bound (|28)) becomes successively weaker 
as the number of terminals is increased. In particular for 
n ^ oc, Pmax(^) is cqual to Pq, whenever \x\ > 1 . 


C. Bounding power by efficiency 

It is instructive to express the results obtained so far 
also as a function of the efficiency defined in (| 8 |) . To this 
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FIG. 5. Bounds on the power -Pmax (? 7 , x) of the n-terminal 
model as a thermoelectric heat engine in units of Pq and as 
functions of the asymmetry parameter x and the efficiency rj 
for n = 3,4 and n ^ oo. The white regions in the plots are 
forbidden by the bound (fT5l) . 


end, we fix the normalized efficiency 77 = 77/770 by putting 


{jl^qp ffi ^pq) ^ {V^qp "b ^pq)‘^ ^V^pp^qq 

‘^Lpp 


(30) 

Inserting this 7 into the general expression (| 9 |) for the 
power, replacing Lqq with its upper bound dm and ex¬ 
pressing the result in terms of the parameters x and y 
yields a bound whose rather involved ex¬ 

pression is given in the appendix. By maximizing this 
function with respect to y, we obtain the upper bound 
Pmax(^, 7) on power for given efficiency 77, asymmetry pa¬ 
rameter X and number of terminals 77., which is plotted in 
Fig. [5l We find that Anax(7,^) vanishes linearly when¬ 


ever 77 approaches its previously identified upper bound 
m- Moreover, in the limit 72 ^ 00, for which the bound 
becomes the weakest, we end up with the re¬ 
markably simple expression 


PmayiijliP) — 4 Po 


r?(l - 7 ]) 
fj — rf Ix^ 


for \x\ > 1 
for |x| < 1 ’ 


(31) 


which becomes independent of x for |x| > 1 . In par¬ 
ticular, m shows that the power must vanish at least 
as ( 4 Po/ 77 c)( 7 c ~ v) whenever 77 approaches the Carnot 
value. We can therefore conclude that the new constraint 
dm rules out the option of finite power at Carnot effi¬ 
ciency for any system that can be described by a multi¬ 
terminal model with an arbitrary number of probe ter¬ 
minals. This insight is our third main result. 


D. Power at maximum efficiency 


After having determined the maximum power of the 
Ti-terminal model as a heat engine for given efficiency 77, 
as another quantity of particular interest, we will now 
investigate power at maximum efficiency. To this end, 
we recall that the maximum efficiency (pTj) of a thermo¬ 
electric generator in the linear response regime is found 
by optimizing the general expression (0 with respect to 
7 under the condition P > 0 . Inserting the thus deter¬ 
mined 7 into ([9j) gives the power at maximum efficiency 


P*{x,y) 


X y/l + ^ - 1 

a/i + y + y +1 


(32) 


By estimating Lqq in terms of its upper bound ([T6|) and 
eliminating y in favor of r]max using m we obtain the 
bound 


P (x, 77max) 


4Po7max 


/hn{x^(^X 7max) 4(r77maxA 
V hn{x){x^ - ) 

(33) 


on the power output at maximum efficiency, where 

7max = 7max/7C* 


Fig. [6] shows plots of the function P*(x, 77 niax) for 
71 = 3, 71 = 4 and 71 ^ oc. For any finite 71, this 
function shows a nontrivial global maximum. Specifi¬ 
cally, we find the maximum 0.726Po at {x^fjmax) — 
(1.132,0.422) for 71 = 3 and the maximum 0.779Po at 
(^,7max) — (1.373,0.432) for 71 = 4. For 71 ^ 00, the 
global maximum converges to Pq. However, this value 
can be reached only asymptotically for fjma^ = 1/2 and 
X ^ ±00. Moreover, we find that, independent of the 
number of terminals 71 and the value of the asymmetry 
parameter x, the power at maximum efficiency vanishes 
like ( 4 Po/r?c)(??max(a;) - r?max) when r/^ax saturates its 
upper bound m- 
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FIG. 6. Plots of the bound (|3^ on the power at maximum 
efficiency of a n-terminal thermoelectric generator in units of 
Po and as functions of the asymmetry parameter x and the 
maximum efficiency ?7max for n = 3, n = 4 and n ^ oo. 

IV. CONCLUDING PERSPECTIVES 

Strong numerical evidence let us to identify the new 
constraint ([T 6 |) on the Onsager coefficients describing the 
coupled heat and particle transport in the presence of 
a magnetic field within the paradigmatic multi-terminal 
model for a thermoelectric generator. We emphasize that 


our numerical scheme uses only the sum rules and the 
structure m of the matrix blocks the primary Onsager 
matrix in (|2]) is composed of. In particular, the only prop¬ 
erty of the function f{E) showing up in (|3]) that enters 
our calculations is f{E) > 0 for any E > 0. Therefore, 
although here provided for a quantum system, with a 
properly adjusted normalization constant Nqq the con¬ 
straint ([ 16 ]) and all of its consequences likewise hold for 
classical models like the railway switch mod el and the 
recently proposed classical Nernst engine for which 
f{E) becomes proportional to the Boltzmann factor. 

In the second part, we used the constraint m to de¬ 
rive an upper bound on the power output (|28|) as a func¬ 
tion of the asymmetry parameter x that quantifies the 
extend to which an externally applied magnetic field lo¬ 
cally breaks the time reversal symmetry of the dynam¬ 
ics inside the scattering region. In the symmetric case 
X = 1 and for /i = 0 our bound assumes the value 
Anax(^ = 1) = {l/2A)7r^klAT^/h with AT = Th - T, 
which is about a factor of 1.3 larger than the bound 
— 0.0321V7r^/c^AT^//i for a single conduction 
channel, i.e., V = 1 , which was derived by Whitney 
[mil for a quantum coherent conductor in the non¬ 
linear regime for vanishing magnetic field. The deviation 
can be explained by the fact that our numerical algorithm 
does not exploit all features of the primary Onsager ma¬ 
trix that result from the structure of formula ([3]) . There¬ 
fore, the constraint m and the bounds it implies should 
not be necessarily regarded as tight, even though they 
seem to be quite good as the comparison with Withney’s 
result for a special parameter value shows. 

Finally, we analyzed the consequences of m for the 
relationship between power and efficiency of the multi¬ 
terminal model as a thermoelectric generator. We found 
that, for any n, the maximum power for given 77 must 
vanish linearly when 77 approaches its upper bound (fT5]) . 
Moreover, it turned out that, for 77 , ^ 00 , the constraint 
m implies the bound m, which is independent of the 
asymmetry parameter x for \x\ > 1 and comprises the 
bound Tmax(^ 5 ^) for any finite n. Therefore, the as¬ 
sumption that any type of inelastic scattering or inter¬ 
action between electrons can be mimicked by a sufficient 
number of probe terminals attached to the scattering re¬ 
gion, tempts us to speculate that, up to the normaliza¬ 
tion constant Pq , the bound (|3T]) is universal beyond the 
multi-terminal model. At the current stage, however, 
an algebraic proof of ([T 6 |) within this setup as well as a 
putative derivation of ([3T]) from first principles without 
reference to a specific model class constitute challenging 
topics for future research. 
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Appendix: Calculation of the maximum power for given efficiency 

The function P±{r],x,y) obtained by the procedure described in section IILC after equation (|3Q|) reads 


P±{v^x,y) = - 4 Po 1 - 


y 


yy [x + T] l + / (x + yP x{l ^ y)y 


hn{x) J l^y \ 2x 


y X 


y 


where y = ylyc denotes the normalized efficiency and y is subject to the constraint 

Axy 1 / \ r 7 / \ 4x7/ 

^ y ^ hn[x) tor X > 0, hn{x) < y < 


(x — y) 


(x — yY 


for X < 0. 


(A. 1 ) 


(A.2) 


Since it is readily seen that P+(7/,x,7/) > P_(7/,x,7/) and we are interested in an upper bound on power, from here 
on, we will only consider P+(?7,x,7/). With respect to y, this function assumes its maximum at 


, _2x + {2x - ri)hn{x) + {2x + {x - ri)hn{x)) 

y [x,r]) = T] - 2 — o - I - \2rT~\ -• 

x^ — Zxy -\- [y — x)^rin[x) 

Consequently, the function Pmax(^,^) discussed in section IILC and plotted in Fig. [5] is given by 

Pmax{V,x) = P+{7],X,y*{x,7])). 

We note that, in the limit ri ^ oc, (Q reduces to 

4X7/ r I I ^ 1 47/ 


(x - 1)(1 + X - 27/) 

and (jA. 4 p assumes the simple form (jST]). 


for |x| > 1, 


X — X^ — 27/ + 2x7/ 


for |x| < 1 


(A. 3 ) 


(A. 4 ) 


(A. 5 ) 
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